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Abstract 

We present recent results for heavy-flavor observables in nucleus-nucleus collisions at LHC energies, obtained with 
the POWLANG transport setup. The initial creation of cc and bb pairs is simulated with a perturbative QCD ap¬ 
proach (POWHEG-i-PYTHIA); their propagation in the medium (created in the nucleus-nucleus or in proton-nucleus 
collision) is studied with the relativistic Langevin equation, here solved using weak-coupling transport coefficients. 
Successively, the heavy quarks hadronize in the medium. We compute the nuclear modification factor and the elliptic 
flow parameter of the final D mesons both in nucleus-nucleus and in (for the first time, in the POWLANG setup) 
proton-nucleus collisions and compare our results to experimental data. 
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1. Introduction 

The purpose of our work is to provide a comprehen¬ 
sive setup for the study of heavy-flavor observables in 
high-energy hadronic and nuclear collisions, from the 
QQ production in hard partonic processes to the detec¬ 
tion in the experimental apparatus. 

In a series of papers llll|2l|3]ll over the last few 
years we developed a complete setup (referred to as 
POWLANG) for the study of heavy flavour observ¬ 
ables in high-energy nucleus-nucleus (AA) collisions, 
describing the initial hard production of the QQ pairs 
and the corresponding parton-shower stage through the 
POWHEG-BOX package fSjIfll and addressing the suc¬ 
cessive evolution in the plasma through the relativistic 
Langevin equation. Here, following Ref ID, we sup¬ 
plement our numerical tool by modeling the hadroniza- 
tion of the heavy quarks accounting for the presence 
of a surrounding medium made of light thermal par- 
tons feeling the collective flow of the local fluid cell. 
Moreover, we present our first (preliminary) results for 
proton-nucleus (pA) collisions. 
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2. Heavy flavour in proton-proton collisions 

Because of their large mass, the initial production of 
c and b quarks is a short-distance process described by 
perturbative QCD (pQCD). Eor this purpose we rely on 
a standard pQCD public tool, namely POWHEG-BOX, 
in which the hard QQ event is interfaced with a shower 
stage described by PYTHIA IT), to include the effects 
of Initial- and Pinal-State Radiation ||Tl|2l. 

Experimental data obtained in pp collisions can be 
exploited to validate the theoretical calculations used 
to simulate the initial hard QQ production. In our 
setup, the heavy quarks are created in pairs by the 
POWHEG-BOX event generator. Eventually, heavy- 
quark hadronization and the final decays of the D (B) 
mesons are simulated with PYTHIA, which is also used 
to describe the parton shower stage. 

In Fig. [T] we show our results for D-h azimuthal cor¬ 
relations compared to preliminary ALICE data MB, 
for two different p^-intervals of the charmed meson. In 
our simulation and D* are used as trigger parti¬ 

cles and the fight hadrons are limited to charged pions 
and kaons, protons and antiprotons, excluding the weak 
decays of A and K^. Any D-meson is correlated with 
all the fight hadrons (except its own decay products) 
created in the same event. The near-side peak takes 
contribution both from correlations present at the par- 
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Figure 2: The R^a of mesons in central Au-Au collisions at 
= 200 GeV. POWLANG results obtained with HTL transport 
coefficients and a decoupling temperature Tdec = 170 MeV are plot¬ 
ted. Also shown for comparison is the extreme case of full kinetic 
thermalization of D mesons. Theory curves are compared to STAR 
data fin 


Figure 1: Azimuthal D-h correlations in pp collisions at ^/s = l TeV 
for different pr-cuts compared to preliminary ALICE data (SE). 


tonic level (from QQ pairs arising from gluon splitting) 
and from hadrons coming from the fragmentation of the 
same string of the parent heavy quark. Our results in¬ 
clude also the simulation of the Underlying Event due to 
Multiple Parton Interactions, performed with PYTHIA 
6.4, which gives rise to the pedestal observed in Fig.[T] 

3. Heavy flavour in pA and AA collisions: setups 


flow in the final state. With a Monte Carlo calculation 
we have simulated a large number of initial configura¬ 
tions and we have obtained a realistic estimate of the 
initial average eccentricity. In a given initial configu¬ 
ration, with each collision occurring at the transverse 
coordinate x, (the single nucleons being randomly lo¬ 
cated with a distribution probability given by a realistic 
nuclear density), the entropy density in the transverse 
plane can be defined as 


i(x) = 


K 


^coU 


^exp 

i=l 


(x-x,)^ 

2cr2 


( 1 ) 


In nuclear collisions the cc or bb production given by 
the pQCD calculation must be modified in two ways; /) 
the nuclear parton distribution functions should be cor¬ 
rected for shadowing or antishadowing effects (we have 
adopted here the EPS09 scheme flOl ); ii) the collid¬ 
ing partons acquire, on average, a larger transverse mo¬ 
mentum during the crossing of the two nuclei (Cronin 
effect), which can be described by a Glauber calcula¬ 
tion II 3 . 

The heavy quarks created in nuclear collisions prop¬ 
agate in a strongly interacting and non-static medium, 
whose properties and evolution is described through hy- 
drodynamical calculations, performed with the viscous 
2-1-1 code of Ref. ifTTl (for the AA case), or with ECHO- 
QGP lIT^ (for the pA case, in 2-1-1 dimensions for sim¬ 
plicity). 

While for the initial conditions in AA collisions one 
can adopt the smooth results provided by an optical- 
Glauber calculation, in a pA collision the event-by- 
event fluctuations in the initial state are extremely im¬ 
portant, since they are the main source of anisotropic 


where cr is the “smearing parameter” used to esti¬ 
mate the effective size of the nucleons. The initial ec¬ 
centricity, which translate into a non-vanishing elliptic 
flow lfT3l . can be then evaluated as (the brackets denote 
the average in the transverse plane, with the entropy 
density in Eq.[2as a weight) 

{/ + x 2 j 

Since a full event-by-event simulation with our hy- 
droH-transport setup would be really demanding, for 
the huge computing and storage resources required, we 
evaluate a realistic average background as follows; for 
a given centrality class we average over all the events of 
the considered percentile (based on the number of N^^), 
after rotating each of them by the event-plane angle 1 /^ 2 , 
so that an initial condition representative of the average 
eccentricity 62 can be obtained and meaningful predic¬ 
tions for the elliptic flow can be provided. 

The propagation of the heavy quarks in the plasma 
is modeled as a Brownian motion by employing a rel- 
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Figure 3: The V 2 of D mesons in Pb-Pb collisions at = 2.76 

TeV. POWLANG results (with HTL transport coefficients) with in¬ 
vacuum and in-medium HQ fragmentation at the decoupling temper¬ 
ature Tdec = 155 MeV compared to ALICE data 1171 in the 30-50% 
centrality class and to the limit of kinetic thermalization. 

ativistic Langevin equation, with the transport coeffi¬ 
cients KijTip) (representing the average squared longi¬ 
tudinal/transverse momentum per unit time exchanged 
by the quark with the medium) obtained within a weak- 
coupling scheme (perturbative QCD and Hard Thermal 
Loop approximation CHIU[3) 

Finally, after the propagation through the dense me¬ 
dium, the heavy quarks decouple and hadronize. In the 
pp case, the hadronization is performed with PYTHIA, 
while for the nuclear collisions we introduce a new sim¬ 
ple model to take into account the effect of the thermal- 
ized medium. Once a heavy quark Q, during its stochas¬ 
tic propagation in the fireball, has reached a fluid cell 
below the decoupling temperature Tdec, it is forced to 
hadronize. One extracts then a light antiquark (up, 
down or strange, with relative thermal abundances dic¬ 
tated by the ratio m/Tdec) from a thermal momentum 
distribution corresponding to the temperature Tdec in the 
Local Rest Frame (LRF) of the fluid; the local fluid four- 
velocity allows one to boost the momentum of 
from the LRP to the laboratory frame. A string is then 
constructed joining the endpoints given by Q and 
and is then given to PYTHIA 6.4 to simulate its frag¬ 
mentation into hadrons (and their final decays). In case 
the invariant mass of the string is not large enough to al¬ 
low its decay into at least a pair of hadrons the event is 
resampled, extracting a new thermal parton to associate 
to the heavy quark. 

We would like to stress that this procedure to simu¬ 
late in-medium fragmentation is new, since in our previ¬ 
ous works the same mechanism as in pp collisions (i.e. 
in-vacuum fragmentation) was employed. With this im¬ 
provement we can first of all provide a realistic estimate 
of the role of the thermal light quarks to explain peculiar 


features of the D meson spectra at low and moderate pj', 
secondly, the complete information on all the final state 
particles arising from the fragmentation of the strings 
allows us to provide theory predictions for observables 
like D-h, e-h, ... correlations to be compared to 

existing data and possibly used as a guidance to future 
measurements. 

Notice that further possible interactions in the hadro¬ 
nic phase, which might enhance the elliptic flow, are 
here neglected; we plan to include them in our future 
work. 

4. Results in pA and AA collisions 

In this section we present a selection of our results, 
compared with experimental data obtained at RHIC and 
LHC energies. 

In Fig.l^we show some POWLANG outcomes for the 
Raa of D^mesons in central (0- 10%) Au-Au collisions 
at y/SffN = 200 GeV. HTL transport coefficients are em¬ 
ployed and the difference between the two hadroniza¬ 
tion schemes (here taken to occur at Tdec = 170 MeVQ, 
either with vacuum fragmentation or with in-medium 
string fragmentation, are clearly visible; in the second 
case the radial flow of the light thermal parton leads to 
the development of a bump around pr ~ 1.5 GeV in 
qualitative agreement with the experimental data. Also 
shown for comparison is the result for the limiting sce¬ 
nario in which charmed particles reach full kinetic equi¬ 
librium with the medium ||4] . Our results are compared 
to STAR datalH. 

In Fig. we address the V 2 of D mesons in semicen¬ 
tral (30-50%) Pb-Pb collisions at - 2.76 TeV at 
the LHC. The effect of the new procedure for in-medium 
hadronization through string fragmentation is clearly 
visible; while POWLANG outcome with standard in¬ 
vacuum fragmentation of charm largely underpredicts 
the observed V 2 , the additional flow acquired from the 
light thermal partons move the theory curves with in¬ 
medium hadronization closer to the experimental data 
measured by the ALICE Collaboration IfTTl . The value 
chosen for the decoupling temperature is Tdec = 155 
MeV; we observed that V 2 is more sensitive to Tdec than 
other observables, the lower value seeming to be pre¬ 
ferred by the data; this agrees with the expectation that 
the elliptic flow needs more time to fully develop with 
respect to the quenching of the spectra. This plot also 


*We notice that the nuclear modihcation factor Raa is not very 
sensitive to the value of Tdec^ with a smaller value (155 MeV) the 
results ai‘e quite similar (5) 
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Figure 4: The V 2 of D mesons in p-Pb collisions at = 5.02 TeV. 
POWLANG results (with HTL transport coefficients) with in-medium 
HQ fragmentation, for cr = 0.2 and 0.4 fm. The V 2 coefficients for c- 
quarks alone are plotted for comparison. 



Figure 5: The nuclear modification factor of D mesons in p-Pb col¬ 
lisions at = 5.02 TeV. POWLANG results (with HTL transport 
coefficients) with in-medium HQ fragmentation, for cr = 0.2 and 0.4 
fm. Experimental data from Ref. GS 

shows how a full kinetic thermalization up to large val¬ 
ues of pt is disfavoured by the data. 

Finally, we show some results for pA collisions. In 
Fig. 12 we show the elliptic flow parameter V 2 for D 
mesons in p-Pb collisions at yfs^ - 5.02 TeV. We 
have plotted the results for two values of the smearing 
parameter cr. The dotted lines show the V 2 for c-quarks 
alone; it is evident the importance of the contribution 
of the light quarks to the collective flow of the charmed 
mesons. In Fig. we plot the nuclear modification fac¬ 
tor for D mesons: we note that the POWLANG setup 
predicts a peak at pj- ~ 3 h- 4 GeV, with a systematic un¬ 
certainty arising from the value of the initial smearing 
parameter. 

We plan to extend this study by considering trasport 
coefficients evaluated on the lattice and to apply the 
same setup also to d-Au collisions, experimentally mea¬ 
sured at RHIC. 


5. CONCLUSIONS 

The simple model to describe heavy quark hadroniza- 
tion in the presence of a hot deconhned medium (a 
Quark-Gluon Plasma) has considerably improved the 
agreement of the POWLANG results with the experi¬ 
mental data at RHIC and LHC energies. In particular, 
results for the nuclear modihcation factors and the el¬ 
liptic flow of D mesons, for pA and AA collisions, have 
been presented. 

With this setup, it is also possible to study heavy- 
flavour correlations in heavy-ion collisions, allowing to 
get additional information on the heavy quark interac¬ 
tion with the medium. 
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